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Abstract 

Flows on the moduli space of the algebraic Riemann surfaces, preserving the periods 
of the corresponding solutions of the soliton equations are studied. We show that 
these flows are gradient with respect to some indefinite symmetric flat metric arising 
in the Hamiltonian theory of the Whitham equations. The functions generating these 
flows are conserved quantities for all the equations simultaneously. We show that for 
1+1 systems these flows can be imbedded in a larger system of ordinary nonlinear 
differential equations with a rational right-hand side. Finally these flows are used 
to give a complete description of the moduli space of algebraic Riemann surfaces 
corresponding to periodic solutions of the nonlinear Schrodinger equation. 
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Introduction. 



In the periodic theory of sohton equations we often meet the following problem. 
The direct scattering transform is well-defined for purely periodic potentials. 
In the quasiperiodic case the spectrum may have a much more complicated 
structure (for example it may be similar to the Cantor set [3]) and in general 
it is impossible to associate a Riemann surface to a quasiperiodic potential. 
When we consider the inverse problem, general algebraic Riemann surfaces 
generate quasiperiodic solutions. 

Periods of solutions corresponding to a given algebraic Riemann surface can 
be expressed in terms of some Abelian integrals. Thus the Riemann surfaces 
corresponding to periodic solutions with a given period form a transcendental 
submanifold in the moduli space of Riemann surfaces of given genus. (More 
precisely we have to consider the moduli space of Riemann surfaces of given 
genus with a set of marked points and jets of local parameters in these points.) 

In the present article we use the following approach to study these submani- 
folds. It has been shown by one of the authors [21] that to any period preserv- 
ing deformation corresponds a meromorphic differential on the spectral curve 
with rather strong restrictions on the singularities. Each such differential gen- 
erates a flow on the moduli space of Riemann surfaces leaving invariant the 
isoperiodic submanifold. Earlier such flows were introduced by Krichever [12] 
for more general meromorphic differentials in connection with the topologi- 
cal quantum field theory and Whitham equations. In our paper the period 
preserving flows are studied more detailed in connection with the problem of 
characterization of the periodic solutions, especially for the 1 + 1 systems like 
KdV and NLS. Flows analogous to those studied in the present paper were 
introduced also by Ercolani, Forest, McLaughhn and Shina [6] in the theory 
of nonintegrable periodic perturbations of the soliton equations. 

These flows can be naturally extended to general finite-gap solutions. In this 
case they preserve the group of frequences. 

In the 1-1-1 case these fiows can be written as ordinary differential equations 
for the branch points but the right-hand side of of these equations contains 
Abelian integrals. We show that these flows can be essentially simplifled by 
extending the conflguration space. If we add to the branch points the zeroes of 
the quasimomentum differential as additional variables we obtain a system of 
ordinary differential equations with a rational right-hand side. These fiows are 
much more convenient both for analytic studies and numerical simulations. 
Differential equations on the zeroes of the quasimomentum differential were 
written earlier in [12]. But our observation is that by combining these two sets 
of parameters together we get differential equations with a rational right-hand 
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side which are much simpher. 

In some sense this system of ODE's is completely integrable, i.e. it can be lin- 
earized in terms of some Abelian integrals. But it does not mean integrability 
in the Painleve sense because the structure of the branch points may be rather 
complicated. 

It is possible to introduce local coordinates on the moduli space given as 
the values of the quasimomentum function in the stationary points of this 
function. These coordinates are usually multivalued. Such coordinates were 
also used by Krichever [12], [13], [15]. In some specific cases like the real 
KdV and the defocusing nonlinear Schrodinger equation they are single- valued 
global coordinates. For the real KdV case similar coordinates were used by 
Marchenko [16]. In this case the moduli space was also investigated by McKean 
and Van Moerbeke [17] with the help of a different coordinate system. 

We show that if we change these coordinates the variations of the correspond- 
ing Riemann surface are also described by ODE's with a rational right-hand 
side. On the basis of these equations a numerical program for studying isope- 
riodic deformations has been written^. 

Isoperiodic flows are rather naturally connected with the Whitham equations 
in the Flaschka- Forest-McLaughlin form [8]. The Hamiltonian theory of the 
Whitham equations was constructed in some specific cases by Forest and 
McLaughlin [9] in terms of action-angle variables and for general situation 
by Dubrovin, Novikov [5] in the local differential-geometrical form and later 
developed by scientists from the Novikov's group. Novikov pointed out to the 
authors that these flows should be gradient in the same fiat Riemann metric 
which arose in the Hamiltonian theory of the Whitham equations. Here we 
give the proof of this conjecture. A flow is called gradient if its vector field 
is defined as the gradient of a function with respect to some symmetric Rie- 
mann metric. Let us recall that Hamiltonian flows are defined as gradients of 
Hamiltonian functions in a skew-symmetric metric. 

In our case it is important that this symmetric metric is not positive definite 
and even degenerate in some points. If we consider a gradient flow in a positive 
definite metric it never conserves the function generating this flow. But in 
our case the differentials of the functions generating isoperiodic flows span a 
maximal isotropic subspace of our metric. Thus all these functions are integrals 
of motion and our gradient system is similar to an integrable Hamiltonian 
system. In fact there also exists a Hamiltonian formulation of isoperiodic flows 
but the corresponding skew-symmetric form is nonlocal with respect to the 
branch points. 



^ The idea of the program interface was suggested to the authors by A. Bobenko. 
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We use this technique in the third part of our paper to study the moduh space 
of all algebraic Riemann surfaces of finite genus generating periodic sohitions. 
We fix attention to the nonhnear Schrodinger equation and show that there 
exists one-to-one correspondence between such Riemann surfaces and some 
graphs together with the flow parameters mentioned above. The multivalued- 
ness of these parameters corresponds to different topologies of these graphs. 
Using these graphs it is convenient to describe the branch points of these coor- 
dinates. To formulate the hypotheses about the structure of the moduli space 
proved later we used some numerical simulations. 



1 Flows preserving the periods. 

We shall study some nonisospectral flows on the space of the flnite-gap solu- 
tions of the soliton equations preserving the periods. 

Let us recall some basic constructions from the periodic sohton theory [23], 
[18]. For deflniteness let us consider the periodic Korteweg-de Vries equation. 

The L operator for the Lax representation for KdV is the one-dimensional 
Scrodinger operator 



If u{x) is periodic with a period P (without loss of generality we shall assume 
P — 1 such that u{x + 1) = u{x)), then L commutes with the shift operator 
'S'a; {Sxf){x) — f{x + 1). The common eigenfunction ^(A, x) is called Bloch- 
Floquet function 



For generic A we have two eigenfunctions of with the eigenvalues /x+ and 
In fact /i+(A) and /i-(A) are two branches of one function holomorphic on 
a hyperelliptic Riemann surface Y, which is called the spectral curve. Finite- 
gap potentials [19] correspond to algebraic Riemann surfaces (i.e. surfaces 
with flnite number of branch points Aq, Ai, . . . , X2g, oo). The function p(A) = 
;^^7= ln(//) is called quasimomentum. It is multivalued. The differential of the 
quasimomentum dp is uniquely deflned by the following properties: 

1) It is meromorphic on Y with only one pole of second order in the point oo 



L = - 



dl+u{x). 



= A*, ^(A,a;+ 1) 



li^{X,x). 
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2) All the periods of dp are multiples of 2n 

^ dp = 2Tlkn, kn ^Ij 

Cn 

over any closed cycle c„. 

In the inverse scattering problem we can consider the quasiperiodic potentials 
as well as the periodic ones. We can start from a hyperelliptic Riemann surface 
with arbitrary finite branch points Aq, Ai, . . . , A2g, (oo is a branch point 
too) . The potentials corresponding to this Riemann surface are parameterized 
by divisors on this surface of degree —g (or equivalently by collections of g 
points 7i, . . . , 7g in y). In general these solutions are quasiperiodic. The basic 
frequences are equal to the periods of the differential dp, uniquely defined by 
the property 1) and 

2') All the periods of dp are real: 

j> dp = 27rfcn, A;^ e M 

over any close cycle c„. 

The potential is periodic with period if and only if for all n 

Pink^ G Z. 

Thus we have the following characterization of the Riemann surfaces corre- 
sponding to periodic potentials (the periodicity properties do not depend on 
the choice of the divisor) : 

There exists a differential dp with the properties 1) and 2). 

We shall study deformations of Riemann surfaces preserving the periods of 
dp. These deformations change the Riemann surfaces. Therefore we call them 
nonisospectral, but they preserve all the x-frequences of the solutions. In par- 
ticular they leave invariant the space of potentials periodic in x with a fixed 
period. Thus we call them isoperiodic. 

For the real KdV condition 2') is equivalent to the to the following: 
2") All the a-periods of dp are equal to zero 

j) dp = Q, 
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where ai . . . . , ag are the basic a-cycles, generated by real ovals. We can extend 
this definition to the complex case and study deformations preserving the b- 
periods of differential dp with normalization 2") where ai, . . . , ag are a-cycles 
from some fixed basis. In contrast to 2') in the complex case this definition 
depends on the choice of the basis of cycles. 

Let us consider the general situation. We assume that the following data are 
given: 

(i) a nonsingular compact Riemann surface Y, 

(ii) a finite number of points yi, . . . ,yn on this Riemann surface and 

(iii) local parameters Wi = 1/ Zi . . . ,Wn = 1/zn near these points. 

Without loss of generality it is convenient to assume that Wi is equal to zero 
at the point y^. 

Now let H be the group of all non- vanishing holomorphic functions / defined 

on some set of the form U\ {oo}, where U C CP"*^ is some neighbourhood of oo 
(which may depend on /). The group multiplication is given by multiplication 
of functions. The Lie algebra of this group contains the vector space h of all 
meromorphic functions on CP^, which are holomorphic on C C CP^. With the 
help of these data (Y, yi, . . . , Zi, . . . , we define a group homomorphism 

L : H"^ Picard group of Y: 

For this purpose we choose some covering 

y = y\{yi,...,y„}UC/iU...UC/„ 

with some small disjoint neighbourhoods of yi,...yn- Then 

L(/i, . . . , fn) is defined with respect to this covering by the cocycle, which 
is equal to z,i i-^ f{z-i) on f/j \ {y^} for all i = 1, . . . ,n. In particular each 
element (/i, . . . , /„) of h" defines a fiow on the Picard group of Y ^ which is 
given by multiplication with L{exp{tfi), . . . , exp(t/„)) for all t e C. These are 
the isospectral flows of the dynamic systems, i.e. these flows do not change 
the Riemann surface. In the sohton theory we have two fixed elements of 
generating the spatial shifts and the fiows generated by other elements are 
associated with generalized time shifts. If one or both of the spatial flows are 
trivial we have a system with a reduced spatial dimension. It is convenient to 
impose conditions on this flows, and therefore to restrict the space of admissi- 
ble data. In doing so the dynamic system may become a Hamiltonian system. 
We will consider three cases: 

finite dimensional case: The both fixed elements of h" induce trivial fiows 
on the Picard group of Y. In this case we have finite-dimensional ODE's and 
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we do not need to discuss periodicity properties. In this case the 'isoperiodic 
deformations' simply preserve the structure of these equations. 

simple periodic case: One fixed element of h" induces a trivial flow and an- 
other fixed element of h" induces a flow with period 1 on the Picard group 
of Y. In this case we have systems with one spatial variable hke the KdV or 
the nonlinear Schrodinger equation and the isoperiodic deformations pre- 
serve the periodicity in this variable. The meromorphic function generating 
the trivial fiow maps our Riemann surface Y to the complex plane and it is 
natural to treat F as a ramified covering of CP^ . The positions of the branch 
points gives us the natural coordinates on the moduli space. 

double periodic case: The both fixed elements of h" induce fiows with pe- 
riod 1 on the Picard group of Y. We have systems with two spatial variables 
like the Kadomtsev-Petviashvili equations and the isoperiodic deformations 
preserve the periodicity in both spatial variables. 

We assume that at all points yi, . . . , y„ at least one of the two fixed elements 
of h" has a singularity. Otherwise we may neglect the corresponding point 
and diminish the number n in our data. It is quite obvious that the fiow 
corresponding to some element / of is trivial, if and only if there exists 
some meromorphic function on Y, which solves the Mittag Leffler distribution 
defined by the cocycle corresponding to /. This means that there exists a 
function holomorphic on F \ {yi, . . . y„} with prescribed singularities fi{zi) + 
0(1) in the points yi, i = 1, . . . n. Furthermore, a fiow corresponding to some 
element / of h" has period 1 if and only if there exists some non-vanishing 
holomorphic function on Y\{yi, . . . y„}, such that the quotient of this function 
divided by the function Zi i— exp{fi{zi)) extends to a holomorphic function 
near yi for alH = 1, . . . n. Hence in the three cases mentioned above we have 
two holomorphic functions A and on F \ {yi, . . . yn}: 

In the finite dimensional Ccise both functions are meromorphic with fixed 
Mittag Lefiler distributions near yi, . . . , 

In the simple periodic case A is a meromorphic function with fixed Mittag 
Leffler distribution near yi, . . . ,yn and is a non- vanishing holomorphic 
function on y \ {yi, . . . ,yn}, such that ln(//) has a meromorphic branch 
with fixed Mittag Leffler distribution near yi, . . . ,yn. 

In the double periodic case both functions A and /i are non-vanishing 
holomorphic functions on F \ , such that ln(A) and ln(/x) 

have meromorphic branches with flxed Mittag Leffler distributions near 

yi,---,yn- 

Hence in all cases these two functions obey some relation of the form 

i?(A,//) = 0, 

where R{-, •) is some holomorphic function on C \ {0} x C \ {0}. In the finite 
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dimensional case R{-,-) may be chosen to be a polynomial with respect to 
A and n and in the simple periodic case only with respect to /i. In fact, in 
both cases the meromorphic function A induces a finite covering of the Rie- 
mann surface. Then there exists an unique polynomial in /i with coefficients 
depending meromorphically on A, which describes the relation between A and 
fi (see [10]). In the double periodic case we will explain later how to define 
this function R{-, ■). This relation determines the Riemann surface Y. In fact, 
the transcendental curve Y defined by the equation -R(A, /i) = may differ 
from the Riemann surface Y only by some singular it iesQ. If n is greater than 
2, it may be convenient to modify our construction and identify the points 
1/1, ...,?/„ to one multiple point. The corresponding singular Riemann surface 
is denoted by An analogous construction defines a group homomorphism 
from H"' into the Picard group of the singular Riemann surface Y'. 

Let us now give our main examples: 

(1) The completely integrable systems, investigated in the papers [1] and [20] 
are examples for the finite dimensional case. 

(2) The KdV equation with periodic boundary conditions is obtained in the 
simple periodic case with n = 1. The fixed elements are f{z) = z and 
f{z) = z"^, z^ generates the trivial fiow and z generates the periodic 
spatial one (see e.g. [18]). 

(3) The non-linear Schrodinger equation with periodic boundary conditions is 
obtained in the modified simple periodic case with n = 2. The fixed fiows 
are (^1,^2) ^ {^1,^2) (the trivial one) and (-21,-22) ^ ("\/— T-21, — v^— T-22) 
(the periodic one) (see [11] and [21]). 

(4) The KP equation with periodic boundary conditions with respect to both 
space variables is obtained in the double periodic case with n = 1 and 
the two functions corresponding to periodic fiows are given hj z z and 
z I— i> az"^ (see [13]). 

The subject of this article is to consider the space of all data, which obey the 
conditions introduced in the three cases above. The methods of this paper 
can be naturally generalized to infinite genus Riemann surfaces with appro- 
priate analytic properties (see e.g. [21]). Two data (Y, yi, . . . , ?/„, zi, . . . , z^) 
and (y, yi, . . . , jjn, zi, . . . , Zn) are called equivalent, if they may be mapped bi- 
holomorphically onto each other. It is quite obvious, that two equivalent data 
either both obey the above mentioned conditions or none of them. In fact, we 
are only interested in the set of equivalence classes of data, which obey one of 
the above mentioned conditions. 



^ We assume, that at all points yi, ■ ■ ■ ,yn the greatest common divisor of the degrees 
of the two fixed functions in h'^ is equal to 1. Otherwise the Riemann surface Y 
may be a non-trivial covering over the curve defined by the equation R{X, fi) = 0. 
^ The definition of the Picard group of a singular Riemann surface is given in [22]. 
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1.1 The finite dimensional case 



Let us now consider the finite dimensional case. Then those data, which obey 
the conditions introduced above are in one to one correspondence with plane 
algebraic curves defined by the equation R{\, fi) — 0. The highest powers 
of this polynomial with respect to A and fi are determined by the two fixed 
elements of h". It turns out that the polynomials R{-,-), which correspond 
to some data obeying the above condition of the finite dimensional case, are 
exactly of the form 



Here the sum is taken over some finite convex set of Nq x Nq. In fact let us 
assume that R{X, n, t) is a onedimensional family of polynomials corresponding 
to some data obeying these conditions. Then the expression 



is a regular meromorphic 1-form of the (in some cases singular) algebraic curve 
defined by the equation R(X,iJ,,t) = 0. This formula shows that this form is 
regular on the domain, where A and fi are finite. The order of the pole of our 
form in the point yj is not greater than the minimal of the orders of A and /j, 
in this point. We shall call such forms weakly singular. If we want to compare 
the different Riemann surfaces corresponding to different values of we may 
choose cither A not to depend on t, such that fi becomes a multivalued function 
depending on A and t or we choose /i not to depend on t, such that A becomes 
a multivalued function depending on /j, and t. In the first case ^dX and in the 
second case ^diJ, is a weakly singular 1-form of our singular Riemann surface. 
This allows us to identify the tangent space of all equivalence classes obeying 
the condition in the modified finite dimensional case with the space of weakly 
singular 1-forms of the corresponding singular Riemann surface. In fact, the 
parameters give a global parameterization of the set of all data obeying the 
conditions of the modified finite dimensional case. The corresponding weakly 
singular 1-forms are given by 



Ro{X, IJ,) + Y1 n^mX^iT, with n^rn £ C. 



m,l 




dR{X,i^,t) (dR{X,fi,t) 
di 1^ dfi 




dX 



djj, — 



d/j, 



dX^ 



dri. 



,m 



dri. 



,m 
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= -AV™ (^^i^) " = AV™ iP^) " dX 

Example 1 Let n be equal to 2 and the two fixed elements of he given by 
{zi,Z2) I— > {zi, Z2) and {zi, Z2) 1— > {zf, —Z2), respectively. 

Then the corresponding polynomials are exactly of the form 

d-l 

R{X, 11) = ^^- X^"^ - ^ riX\ with ri e C. 

1=0 



The weakly singular 1-forms of the plane curves defined by the equation 
R{X, IJ>) — are given by 

5A , c^/^ „ A^ „ 7 , ^ , 
Tr—dix — TT—dX — ——dX with I — 0, . . . ,d — 1. 
ori dri 2/1 



1.2 The simple periodic case 



Now wc consider the modified simple periodic case. It corresponds to solu- 
tions of 1 + 1 systems, which are periodic in one spatial variable. Due to this 
condition the function •) is of the form 



M-l 
m=0 



with some entire functions ro(-), . . . , rM-i(-) of finite order which have essential 
singularities at A = 00. Hence the transcendental curve defined by the equation 
it!(A, //) = has infinite algebraic genus. Let us again assume that i?(A, t) 
is a onedimensional family of entire functions corresponding to some data 
obeying these conditions. The quasimomentum p is the multivalued function 
defined by y=ln(/i). Then the expression 

dp,. dX 

_ dR{\ii,t) f dR{\i,,t) \-\ _ ^r-^ dRiX,i,,t) / dRiX^y'^^ 
dt \ dX J ^ dt \ dfj. J 
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is a regular meromorphic 1-form on Y. If the deformation of it! corresponds 
to a continuous deformation of Y then this differential is also meromorphic 
on Y{t). This formula shows again that this form is regular on the domain 
where A and are finite. The order of the pole of our form in the point yj is 
not greater than the minimal of the orders of A and p in this point. We shall 
call such forms weakly singular. Thus any weakly meromorphic differential 
cu — udX of Y generates an isoperiodic flow defined by 



(This representation is rather similar to the form in which the Whitham equa- 
tions for the KP equation were written [14]). Let A^ be a simple branch point 
of our covering Y — > CP^. Then from (1) it follows: 



(Similar representation is widely used in the Whitham theory [8]). If all the 
branch point are simple we get a system of ODE's but with a rather compli- 
cated right-hand side. If we want to compare the different Riemann surfaces 
corresponding to different values of t, we may choose either A not to depend 
on such that becomes a multivalued function depending on A and or 
we choose not to depend on such that A becomes a multivalued function 
depending on /i and t. In the first case and in the second case ^dp is a 
weakly singular 1-form of our singular Riemann surface. Thus we arrive at the 
following conclusion (see [21]): The tangent space to all data corresponding to 
solutions which are periodic in one variable with given period is in one-to-one 
correspondence with the space of weakly singular 1-forms of the singular Rie- 
mann surface Y{t). If we also admit Riemann surfaces of infinite genus, the set 
of all functions i?(-, ■), which correspond to data obeying the condition of the 
simple periodic case, is again an affine space. But then we have to take care 
of the boundary conditions near the points yi, . . . , y„. The normalization Y of 
our transcendental curve Y (the surface obtained by removing all the singu- 
larities) is now algebraic and coincides with the algebraic curve which plays 
the role of the spectral data in the finite-gap theory. Consider the following 
meromorphic function on Y 




(1) 



dt 



^(Afc) 
dp{Xk) 



(2) 



dp 1 d/j, 
dX II dX 



1 dR{X,fi) dR{X,fi) 

'Ji ax ' ax 



We may describe Y by some equation of the form 



Q{k, A) = 0. 
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Let us remark that this representation differs from the standard one. This 
function Q{-j •) is again uniquely defined as a polynomial in k, with coefficients 
depending meromorphically on A, which is even a meromorphic function on 
CP^ X CP^, if the corresponding Riemann surface has finite genus. Let us now 
assume that some regular 1-form cu is given by 

uj — vdX — vKT^dp 



with some meromorphic function v on our Riemann surface. Then we may 
calculate the derivative of the function Q{-^ ■) with respect to the corresponding 
flow-parameter t. In order to do this we choose A to depend not on t, such 
that both functions n and ii are multivalued functions depending on A and t. 
By definition we have the equation 

dp 



Equation (1) takes the form: 



d'^p Ok du . , , ddp , 

— = — or equivalently --— = dv. (3j 



dtdX dt dX ^ " dt 



This representation is similar to the Flaschka-Forest-McLaughlin represen- 
tation [8] of the Whitham equations. Equation (3) has been derived in the 
theory of the nonintegrable periodic perturbations of the soliton equations in 
[6]. Similar equations for meromorphic u (which may be not isoperiodic) were 
written in [12] formula (7.68)). Since we have chosen A not to depend on t we 
may calculate with this formula the derivative of Q with respect to t: 

dQ{K,X,t) dKdQ{K,X) dvdQ{K,X) 
dt dt dn dX dn 



If we now could find 

(i) some parameterization of a set of functions Q{-^ ■), which contains the 
functions describing the data obeying the conditions in the modified simple 
periodic case and 

(ii) all meromorphic functions such that vdX is a weakly singular 1-form 
of the Riemann surface described by Q{k, X) — 0, 

then this equation defines some fiows on the space of parameters, which leaves 
invariant the subset of those values of the parameters, which correspond to 
data obeying our condition. If furthermore we 
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(iii) find one function Q{-,-), which corresponds to some data obeying our 
conditions, 

then we may determine a whole family of such functions Q{-,-)- 

In the previous considerations u was an arbitrary weakly singular form. We 
now present a natural set of parameters. Consider a canonical basis of cycles 
ai, . . . ,ag,bi, . . . ,bg, such that the intersection number between Oj and 6j is 
equal to 1 for alH = 1, . . . , and vanishes on all other combinations (see [10]). 
Then the first homology group is isomorphic to the free abelian group 

Zoi © . . . ® ZOg e Z6i ® . . . © Zbg. 



Since the intersection form is a non-degenerate antisymmetric Z- valued form, 
we may associate to the function /j, an element m of the first homology group, 
such that the integral of the 1-form dp is equal to 27i times the intersection 
form with m. By changing the basis we may achieve that 

^ (ip = for alH = 1, . . . , 5f m e Zoi © . . . © ZOg. 

ai 



Then there exists a unique basis uji,...,0i)g of holomorphic 1-forms of our 
compact Riemann surface, such that 

ujj = 27rV^5ij for all i, j e {1, . . . , g}. (4) 



Now we can define g fiows on the space of compact Riemann surfaces of genus 
g, which obey the condition of the simple periodic case: 

-J-dX — T^dp = uji for all i = 1, . . . , o. (5) 



Of course the canonical basis extends uniquely to some neighbourhood of 
each of these Riemann surfaces. But globally the space of Riemann surfaces 
together with some canonical basis will be a non-trivial covering over the space 
of Riemann surfaces. It is quite easy to see, that the parameters ti, . . . ,tg 
locally may be chosen to be given by 

ti = j> pdX for all i = 1, . . . , gr. (6) 



These parameters depend only on the choice of the canonical basis obeying 
the foregoing condition. Moreover, this shows that all these fiows commute. 
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Now it is obvious that there is a number which does not change under this 
flows: 

The greatest number m G N, such that m/m is an element of the first homology 
group. We will see in the example of the nonlinear Schrodinger equation that 
this number classifies the connected components of the set of data obeying 
the conditions of the modified periodic case. 

Remark 1.1 These flows can be extended to Riemann surfaces corresponding 
to general finite-gap solutions of soliton equations with one spatial dimension 
(in general they are quasiperiodic) . In this case the Hsoperiodic deformations' 
preserve the group of the frequences of these solutions. So we consider all data 
which obey the following condition in the 

simple quasiperiodic case: one of the fixed elements in generates a 
trivial flow and no conditions on the flow generated by the second fixed 
element are imposed. 

The function A corresponding to the trivial flow maps Y to CP^ so Y is repre- 
sented as a ramified covering of CP^ . The position of the branch points gives 
us local coordinates on the space of such data. 

In this case an analogue of the quasimomentum differential can be defined 
uniquely as an Abel differential dp of the second kind such that: 

diff 1 : dp is holomorphic m y \ (yi U ^2 U . . . U y„) . 

dijf 2: dp — -^^dfk+reg. terms in the point yk where (/i, . . . , /„) is the second 



fixed element of H"^ 
diff 3: Im§^dp = for any closed cycle c. 

n = dp/dX is again a meromorphic function on Y and the foregoing discussion 
of the isoperiodic flows carries over to this more general situation. In fact, the 
right-hand side of (3) is an exact differential so all the periods of the left-hand 
side are equal to 0. Thus all the periods of dp are integrals of motion. 

Let us consider concrete equations in more detail. 

Example 2 ^5 we mentioned above, the KdV equation with periodic boundary 
conditions is an example of the simple periodic case [19]. It is well known (see 
e.g. [17]), that in this case the function R{-, •) is of the form 



R{X,li) = - II A{\) + 1. 
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Hence we may calculate the function n in terms of A and ji: 



^J-ldR{\n) ,dR{\^ji) dA{\) 
i^^— ^^ — /- 



II dX ' djj dX 2/i-A(A)' 

Hence the function Q{-,-) has the form 



// the Riemann surface has finite genus, there are only a finite number of 
zeroes of^^^, where the denominator /^^{X)— A is not zero too. Moreover, by 
definition of the two fixed elements of h: to the function X there corresponds 
the element z ^ and to the function In there corresponds the function 
z ^ z, hence for very large X the function is almost equal to 1/4A. Then 
Q{-, •) is of the form 



Q{k, X) = k — 



2 nf=i(A-ai)^ 



4n?£o(A-A,)' 



where g is the genus, ai, . . . , ccg are the values of the function X at the 2g zeros 
of the differential dp and Aq, . . . , X2g, oo are the values of the function X at the 

zeroes of the differential dX ( these zeroes are of course the branchpoints of the 
covering map induced by X). The standard representation for this Riemann 
surface differs from this one and reads as: 

C^ = 4(A-Ao)...(A-A2,). (7) 



Now the holomorphic 1-forms of the Riemann surface defined by the equation 
Q{k, a) = are of the form 



a; = ^^ (8) 
2^) 



Here o(-) denotes a function Y^fZo otX^ , d{-) denotes the function d{X) — 
Y^=q{X — Xj) and q{-) denotes the function q{X) = nf=i(A — Oii). The flow 
corresponding to this 1-form is given by the differential equation 

2,(A)i?f)-,(A)^f) = (2.(A)^^-o(A)ff)). (9) 
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This flow may also be formulated in terms 0/ A(-); 

dA{X) _ o{X)dA{X) 
dt ~ q{X) dX 



The flow is written as a system of ordinary differential equations (9), hut 
this system is rather complicated, because the coefficients of the polynomial 
q{X) (or, equivalently, the zeroes of q{X)) depend on the branch points in a 
rather complicated way (they can be expressed in terms of some hyperellip- 
tic integrals). But we may consider Xq, . . . , X2g, ao, ■ ■ ■ , ag, as independent 
variables. This means that the differential 

dp = ^^^ dX = KdX. (10) 
2^ 



will have arbitrary periods. Neverthelesson this bigger set of parameters (9) 
still defines an ordinary differential equation. Moreover, the coefficients of 
o(-) may depend on all these parameters in a rather complicated way, which 
is the case in the choice of (4) and (5). So let us present another choice, in 
which (9) simplifies to ODE's with a rational right-hand side. If all the values 
Qfg, Ao, . . . , X2g are pairwise different, we have a basis of holomorphic 
1- forms 

Ui — —- — ^ — dp for all i — 1, . . . , g with some q e C \ {0}. 
X- ai 



In the periodic case the foregoing formula shows, that the values of the function 
A{ai) depend only on the flow parameter corresponding to the form cui. It is 
well know, that A(A) = 2 cos (y/ X — Xoj describes the Riemann surface (of 

genus 0) corresponding to constant potentials. This transcendental curve has 
infinitely many double points at the values ai = (Aq + in)'^ for all i E N. So 
let us choose this curve as the starting point of our flows and normalize the 
1-forms by the condition 

A{ai) = (-1)*2 cosh(ii) for all i e N. (11) 



So we introduce an infinite sequence of flow parameters. But we are interested 
in the case, where only finitely many of these parameters are not equal to 
zero. This is precisely the case of Riemann surfaces of finite genus. In fact, 
generically the genus is equal to the number of flow parameters, which are not 
equal to zero. The corresponding differential equations are given by 

—i^--^ withj^0,...,2g; (12) 

oti Xj — ai 
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with j = 1, . . . , i — 1, i + 1, . . . 



aj — ai 





(13) 



aj — ai 



The normalization condition implies 




nj=o(«i - \j) 



Locally this uniquely defines the coefficients Ci, . . . ,Cg. Due to the normaliza- 
tion condition all these flows commute locally. It is quite easy to see that the 
singularity of these flows at the starting point may be removed. The i-th flow 
opens the i-th gap. For real flow parameters these flows have no singularities, 
but in the complex domain they have very complicated singularities. The for- 
going formula implies that J^jLo Aj — 2 J2i=i o^i does not depend on t. It is well 
known, that this constant may be an arbitrary number inC (it is related to the 
mean value of the periodic potential). The integrals over all generators of the 
first homology group of the form dp are also invariant under these flows. Hence 
we have 2g -\- 1 integrals of motion. Our space of parameters has dimension 
3g -\- 1. This shows that the set of data with genus g obeying the conditions of 
the periodic KdV equation has dimension g + 1. 

Consider now the general finite- gap quasiperiodic KdV solutions. The formulas 
(7), (9), (12) and (13) are valid in this case too. But in the quasiperiodic 
case we have no natural analogues of A{X). Nevertheless, the function p{X) 
is locally well-defined and the normalization condition (11) is solved by the 
local parameters ^/—lp{ai) -\- iir, i — 1, . . . , g of the moduli space of Riemann 
surfaces corresponding to the solutions with a given set of frequencies. Similar 
coordinates were used by Marchenko [16] in the periodic KdV theory and by 
Krichever in the general situation (private communication). All the periods of 
dp are integrals of motion of (9). In the space of all X^, ak, we have a subspace 
defined by: 



and the flows (9) with arbitrary o(A) leave this subspace invariant. Thus it 
is sufficient for us to know one point of this subspace to construct the whole 
family. In fact these conditions are integrals of motion of our equations. 

Remark 1.2 H. Babujian has observed that the right-hand side of (13) with 
Ci = 1 appeared also in the quasiclassical limit of the Bethe ansatz equations, 
which are the Bethe ansatz equations for the Gaudin magnet [2]. The roots of 
this similarity are not clear now. 



Im f dp = 
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Example 3 Let us consider now the non-linear Schrddinger equation as a 
second example of the simple periodic case. It is very similar to the KdV 
equation. So let us only mention the modifications. The function R{-, •) is 
again of the form 

i?(A,At) =//'-//A(A) + l. 

The corresponding Riemann surfaces of finite genus have two covering points 
of \ = oo, which are identified to one multiple point. They are described again 
by the equation 

where g is the algebraic genus of this singular Riemann surface ( this is equal 
to the genus of the normalization plus one), ai, . . . ,ag are the values of the 
function A at the 2g zeros of the differential dp and Ai, . . . , X2g are the values 
of the function A at the zeroes of the differential d\ (these zeroes are of course 
the branchpoints of the covering map induced by X). 

The standard representation for this Riemann surface has the form 

C^ = (A-Ai)...(A-A2,) = d(A). 

Now the regular holomorphic 1-forms of the singular Riemann surface defined 
by the equation Q{k, A) = are of the form 

a; = ^ ,^ dX with Oj e C for all i — 0, . . . , g — 1. 

i=o ^d{X) 

dp = -^LJ=(iA. 
y/d{X) 

If all the values ai, . . . , ag, Ai, . . . , X2g are pairwise different, we have again a 
basis of regular 1-forms 

uji — —- — ^ — dp for all i — 1, ... g with some Ci & C \ {0}. 
A — ctj 

The values of the function A(aj) again only depend on the i-th flow parameter. 
The function A(A) = 2cos(A) describes the Riemann surface (of genus 0) 
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corresponding to the zero potentials. So let us choose this curve as the starting 
point of our flows and normalize the 1- forms by the condition 

A(ai) = (-l)*2cosh(ti) for all i G Z. 



Again we introduce an infinite sequence of flow parameters. But mostly we 
assume all but a finite number of parameters to be equal to zero. This is again 
the case of Riemann surfaces of finite genus. In fact, generically the algebraic 
genus is equal to the number of flow parameters, which are not equal to zero. 
The corresponding differential equations are given by 

dXj Ci . , . ^ 

— = with J = 1, . . . , 2t/; 

oti Xj - ai 

— = — - with J = l,...,i-l,i + l,...g; 

doii 1 > Ci 



The normalization condition implies 

Uj^iiai- ajY' 



Due to the normalization condition all these flows commute locally. It is quite 
easy to see that the singularities of these flows at the starting point may be 
removed. The i-th flow opens the i-th gap. For real flow parameters these 
flows have again no singularities, but in the complex domain they have very 
complicated singularities. The forgoing formula implies that J2jLi Aj— 2 I]f=i 
does not depend on t. Due to the condition that near X = oo p has near A = oo 
two branches of the form ±A + 0(A~-'^) this constant must be equal to zero and 
furthermore, the integral of the 1-form dp along a path connecting the two 
covering points of X = Xq must be of the form ±(2Ao + 2m7r) + 0(Ag ^) near 
the point Aq = oo. Moreover, the integrals over 2g — 2 generators of the first 
homology group of the form dp are also invariant under these flows. Hence we 
have 2g integrals of motion. Our space of parameters has dimension 3g. This 
shows that the set of data with algebraic genus g obeying the conditions of the 
periodic nonlinear Schrddinger equation has dimension g. 
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1.3 The double periodic case 



Let us now consider the double periodic case. Our main example will be the 
KP equation (see [13]). In general the function R{-,-) , which describes the 
relation between the functions A and yU, is not (uniquely) defined. Wc assume, 
that for all A G C \ {0} the sum over the sheets of the covering map induced 
by A of the function converges. In this case we may define the function 
R{X, li) in the following manner: 

R{x, //) = n 

over the sheets of the covering map induced by A 




This choice may be characterized by two properties: 

(i) For all fixed Aq G C \ {0} -R(Ao, ■) is an entire function with genus equal 
to zero. 

(ii) For all A e C \ {0} i?(A, 0) = 1. 

Now let us again assume that R{X, /i, t) is a one dimensional family of such 
functions, corresponding to some data obeying the condition of the double 
periodic case. Then the expression 

— - — din a — ^-—^dlnX — 
dt ^ dt 

dR{X,fi,t) f^^ 9R{X,fi,t) \-' 



dt \ dX 



dR{X,i,,t) f^dRiX^y ^^^^ 



dt \ dji 



is a regular meromorphic 1-form of the (in some cases singular) algebraic curve 
defined by the equation i?(A, /x, t) — 0. This formula shows again that this form 
is regular on the domain, where A and /i are finite. Due to the conditions on 
the functions A and /i this form has poles of order at most one at all the points 
yi,...yn- Hence it is a regular form of the singular Riemann surface, where 
all the points yi, . . . ,yn are identified to one multiple point. If we want to 
compare the different Riemann surfaces corresponding to different values of t, 
we may choose either A not to depend on t, such that /i becomes a multivalued 
function depending on A and t or wc choose fi not to depend on t, such that 
A becomes a multivalued function depending on fi and t. In the first case 
^^dln A and in the second case ^^dln is a regular 1-form of our singular 
Riemann surface. This again allows us to identify the tangent space of all 
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equivalence classes obeying the condition in the modified finite dimensional 
case with the space of regular 1-forms of the corresponding singular Riemann 
surface. Let us now restrict to the case n = 1. In this case the Riemann surface 
Y may be chosen to be the non-singular normalization of the transcendental 
curve described by the equation R{X, /i, t) — 0. Moreover, we restrict ourselves 
to the case in which this normalization is a compact Riemann surface of finite 
genus g. Then there exists a canonical basis of cycles ai, . . . ,ag,bi, . . . ,bg, such 
that the intersection number of Oj and bi is equal to 1 for alH = 1, . . . , (7 and 
vanishes on all other combinations (see [10]). Then the first homology group 
is isomorphic to the free abelian group 

Zai © . . . ® Za^ © Z6i © . . . © Zbg. 

Since the intersection form is a non-degenerate antisymmetric Z- valued form, 
to both functions A and fi we may associate elements 1 and m of the first 
homology group, such that the integral of the 1-forms dlnX and dlnfi is equal 
to 27r\/— 1 times the intersection form with 1 and m respectively. By changing 
the basis we may always achieve that either 

(i) 0? In A = for all i = 1, . . . , g', 

or 

(ii) 4 = for all i = 1, . . . , 5. 

Furthermore, if the intersection number of 1 and m is equal to zero, we may 
even attain, that both conditions are fulfilled. 

Lemma 1.3 If there exists a double periodic non-singular solution corre- 
sponding to the compact Riemann surface Y , then the intersection number 
of 1 and m is zero. 

Proof: Both elements 1 and m can be considered as elements of the Lie algebra 
of the Picard group of Y. Now it is well known that any double periodic 
solution corresponding to the compact Riemann surface Y may be described 
by the family of divisors 

D[s, t) = exp(sl -I- tTa)D with some divisor D of degree g and (s, t) e M^. 

If this solution is non-singular, all divisors of this family are non-special and 
have support inside of y . . . Hence all these divisors may be described 
by an unique gi-tuple of points of Y . Since both fiows are assumed to be 
periodic, the two parameters (s,t) of the family run through (m/z)^. The 
summation over the points of the divisors defines a map 

//,((R/Z)2,Z) ^i7,(F,Z). 
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Under this map the two generators of Hi{{r/z)^ ,Z) are mapped onto 1 and 
m. Moreover, this map respects the intersection number. If the intersection 
number of 1 and m is not zero, the image of the generator of i72((K/z)^, Z) 
cannot be equal to zero. Then the union of the paths of all points of the 
divisors of the whole family covers the whole Riemann surface Y. This is 
a contradiction to the assumption, that all divisors of the whole family are 
non-special and have support inside of F \ {yi, . . . □ 

In the sequel we assume both conditions (i) and (ii) . Then there exists a unique 
basis cui, . . . ,u!g oi holomorphic 1-forms of our compact Riemann surface, such 
that 

j) Uj = Sij for all i, j e {1, . . . , g}. 



Now we can define g flows on the space of compact Riemann surfaces of genus 
g, which obey the condition of the double periodic case: 

^—-^dln a — ^-—^dlnX = uii for alH = 1, . . . , o. 
dti dti ' 



The canonical basis of course extends uniquely to some neighbourhood of 
each of these Riemann surfaces. But globally the space of Riemann surfaces 
together with some canonical basis will be a non-trivial covering over the 
space of Riemann surfaces. It is quite easy to see, that locally the parameters 
ti, . . . ,tg may be chosen to be given by 

ti = ln(/x)(iln(A) for alH = 1, . . . , gi. 

at 



Due to our assumption the multivalued function In/x may be chosen to be 
single valued on each of the cycles ai,...,ag. Hence these parameters are 
locally defined only up to summation of some elements of Z27r\/— 1. If the 
intersection number of 1 and m is zero, these parameters depend only on the 
choice of the canonical basis obeying both conditions (i) and (ii). Now it is 
obvious that three numbers do not change under these fiows: 

(1) the greatest number / e N, such that l/l is an element of the first homology 
group. 

(2) the greatest number m e N, such that m/m is an element of the first 
homology group too, and 

(3) the intersection number between 1 and m. 

Isoperiodic flows can be naturally extended to the general quasiperiodic finite- 
gap solutions. Let y be a compact Riemann surface. Then any holomorphic 
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differential on Y generates a deformation preserving all the x and y frequences. 



2 Gradient flows in the simple periodic case 



In the section 'simple periodic case' the isoperiodic flows were written in terms 
of the branch points of "K Ai, . . . , Ajv as a system of ordinary differential equa- 
tions. These equations are defined for the finite-gap quasiperiodic solutions as 
well as for the periodic ones. In this case these flows are defined on the space of 
nonsingular compact Riemann surfaces such that one of the fixed elements is 
represented by a meromorphic function A and generates a trivial fiow. A maps 
Y to CP^. Hence the zeroes of the quasimomentum arc uniquely defined by 
the branch points and the branch points are the only independent variables. 

The Hamiltonian theory of the Whitham equations [5] is based on the existence 
of the following Riemann metric on our moduli space: 

N 

ds^ = {dXif , with 
1=1 



Qii = res 



A=Ai 



^foralli = l,...7V, ^. _ q for i ^ j. (14) 



(This formula is due to Dubrovin [4]). 

In this section we show that isoperiodic fiows are gradient fiows in the met- 
ric (14). This fact was pointed out to the authors by S.P.Novikov. We shall 
consider the fiows (3) where the differential dp is normalized by: 

rfp = 0, i = l,...,5r. (15) 



Such fiows preserve the 6-periods of dp. 

Theorem 2.1 Let cuj be a holomorphic differential from the standard basis 
(4)- Then the corresponding fl,ow reads as: 

dtj ^ dxr ^ ' 

where 



H.^-l-j>dp. (17) 
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(The metric tensor gki and the inverse matrix g^^ are diagonal g^i — Slgkk, the 
corresponding Riemann metric is flat [5].) 

The isoperiodic property means that the functions Hj are conserved quantities 
for all flows 



The representation (16) is very similar to the standard Hamiltonian represen- 
tation, but we have a symmetric form g''^ instead of a skew-symmetric one in 
the Hamilton theory. 

In this example it is very important that this metric is not positive definite. 
For a positive definite metric the function generating the fiow can not be a 
conservation law. But for the fiows (16) the differentials of the functions dHj, 
j = 1, . . . , (7 lie in the 'light cone' of this metric. Thus all of them are conserved 
quantities for the whole system. From this point of view equation (16) is rather 
similar to integrable Hamiltonian systems. 

Remark 2.2 Equation (16) can he rewritten in a Hamiltonian form with the 
same 'Hamiltonians ' but the corresponding skew-symmetric form has no nat- 
ural representation like (14)- 

Remark 2.3 Introducing the flat coordinates of (14) we explicitly linearize the 
system (16). But the metric g^k has singularities and the flat coordinates have 
rather complicated singularities in these points. The metric g^k is nonsingular 
if and only if dp and dX have no common zeroes. 

Proof of the Theorem 2.1. Following the paper [4] we will introduce flat 
coordinates for the metric (14). 

For this purpose denote the order of the pole of the function A in the point yk 
hj Ik + 1. It is natural to choose the local parameters Wk in the points yk so 
that 



dHj 



— for all j, k. 



dtk 



Consider the following collection of functions on our moduli space: 



Hk = -7^ f dp, tk = 




1 



/ 



pdX, k = l,...,g. 
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Sk = -leSy^^pdX, k = 2,. . .,n. 
dp 



k(wa 



1, . . . , 71, 1, . . . , /q,. 



In the points of general position these functions give us a local coordinate 
system. The metric (14) in these coordinates takes the form 



9 

k=l k=2 



ds^ = {dHkdtk + dtkdHk) + idr^dsk + dskdvk) + 



la 



+ H + 1) 51 dta,kdta,la+l-k- (1^ 



a=l k=l 



From (3) it follows that the flows generated by holomorphic differentials do 
not change the variables r^, s^, ta^^. 

Comparison of (18) and (16) finishes the proof. 



3 The moduli space of the periodic NLS equation 



In this section we want to give a description of the set of data corresponding 
to periodic solutions of the non-linear Schrodinger equation. It is quite easy 
to do the same for the KdV equation. As mentioned above the Riemann sur- 
faces corresponding to such periodic solutions of the non-linear Schrodinger 
equation are described by the equation 

/i^ - fiA{X) + 1=0 



with some entire function A. This function may be described as an infinite- 
sheeted covering A : CP^ — CP"*^ with an essential singularity at infinity. If 
the Riemann surface is of finite genusQ the function In fi is meromorphic on 
both sheets over some open neighbourhood of A = oo. Hence the function 
arccos(A/2) describes all sheets of the covering map A over some neighbour- 
hood of A = oo. Furthermore, this local coordinate near this branchpoint of 
infinite order extends to all but a finite number of the sheets over the complete 
domain A e CP^. Finally for big |A| In /x is of the form ±^/—lX + 0(1/A). 
Hence the global parameter A of the covering space which covering is induced 

^ With the help of the analysis of [21] our discussion may be extended to the infinite 
genus case. 
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by A determined Tip to sign and summation of multiples of 27r\/— 1 by this 
covering map A : CP^ CP^. Later we will label the sheets of this covering by 
the integers and will give a glueing rule for these sheets. The glueing rule will 
uniquely determine the parameter A. This observation proves the following 
lemma: 

Lemma 3.1 Let A i-^ A(A) correspond to a Riemann surface of finite genus of 
a periodic solution of the non-linear Schrddinger equation. Then this function 
is completely determined by all the values of A at the branchpoints of the 
covering map 



In the first section we introduced flow parameters (ii)iez, such that the values 
of A at all branchpoints of the covering map 



are equal to (— 1)*2 cosh(tj). In order to apply the previous Lemma we have to 
add to these values the glueing rule for the covering map. Let us first consider 
the starting point, when A(A) is equal to 2cos(A). In this case we choose the 
following glueing rules: We have infinitely many sheets labeled by some index 
i e Z. For all j G Z the 2j — 1-th sheet and the 2j-th sheet has branchpoints 
at A = —2, oo and the 2j-th sheet and the 2j + 1-th sheet has branchpoints 
at A = 2, oo, and the cuts are chosen as indicated in figure 1. The parameter 
A is uniquely defined by the following choice: On the 2j-th sheet it is equal to 
{2j - l)nV^ for A = -2 and equal to 2j7rV^ for A = 2. 



On the 2j + 1-th sheet it is equal to (2j -|- l)7r-\/^ for A = — 2 and equal to 
2jny/^ for A = 2. 



C ^ C, A 1-^ A(A) 



together with some glueing rule for the infinite number of sheets. 



□ 



cpi ^ cp\ A A(A) 



cut bet'ween the 



cut )etween the 
2j-tli and the 
(2j -- l)-th sheet 



{2j - l)-th and the 
2j-th sh<;et 



A = -2 



A = 2 



Figure 1. 
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With this glueing rule it is quite obvious what happens, if two cuts along a 
common sheet are passing each other. In the following figures we move only 
both cuts along the 2A;-th sheet for some fixed A; e Z and assume that j runs 
through z\{k}. The direction of the movement of the cuts from one figure to 
the next figure is indicated by arrows: 



A = -2 



cut bet'ween the cut between the cut between the cut 

{2j - l)-th and the {2k - l)-th and the 2A:-th and the 2j- 
2j-th sh(!et 2k- :h sheet {2k + l)-th sheet {2j 




Figure 2. 



)etween the 
and the 
l)-th sheet 



til 



A = 2 



cut bet'vMeen 

(2j - 1 
2j-th sh^et 



the cut 
th and the 2k 
{2k 



A = -2 



between the 
;h and the 
+ l)-th sheet 




the cut 
and the 2j 
sheet {2j 



)etween the 
and the 
l)-th sheet 



til 



A = 2 



Figure 3. 



cut bet'wieen 

(2j - 1 
2j-th sh^et 



the cut 
th and the {2k 
{2k 



A = -2 



between the cut between the cut 

l)-th and the {2k - l)-th and the 2j 
+ l)-th sheet 2/c-th shej^t {2j 



A = -1 



jetween the 
and the 
l)-th sheet 



th 

-- 



Figure 4. 
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From the figures 2. to the figure 4. the both cuts along the 2/c-th sheet have 
moved once around each other. Finally the indexes {2k — 1,2k, 2k + 1) are 
permuted to {2k + 1,2k — 1,2k). This shows that the coordinate given by 
the values of A at the branchpoints has a singularity of the form {A2k~i — 
A2fc)^^'^, where = (— 1)^^^"'^2 cosh(t2fc) denotes at the starting point the 
value of A at the branchpoint between the 2k — 1-th and the 2fc-th sheet and 
A2fe = (— 1)^'^2 cosh(t2fc) denotes at the starting point the value of A at the 
branchpoint between the 2A;-th and the 2k + 1-th sheet. This gives a complete 
picture of the branchpoints of the coordinates {ti)i^j^. 

Let us now describe the structure of the covering of the moduli space with 
respect to these coordinates. For reasons of simplicity let us assume that the 
real parts of the values of A at two branchpoints are different, whenever they 
belong to a common sheet. In this case all the cuts may be chosen of the 
form Aj + a/— 1[0, oo]. Hence we have to associate to the values of A at all the 
branchpoints the numbers of the two corresponding sheets. Due to the previous 
Lemma this completely determines the corresponding data. These data may 
be described by a graph with vertices indexed by z G Z corresponding to 
the sheets and links corresponding to the branchpoints and a function from 
the links into the complex numbers, which is equal to the value of A at the 
corresponding branchpoint. These data have to obey the following conditions: 

(i) The graph is connected. 

(ii) All but finitely many vertices are linked exactly with those vertices, whose 
index differs by ±1. 

(ill) The number of links is equal to the number of pairs of the form {i,i + 
1) with z G Z. Due to condition (ii) this has a precise meaning. 

(iv) For large \i\ the number corresponding to the link connecting the i-th 
and the i + 1-th vertex is equal toQ (— 1)*2 

It is quite obvious how this description may be extended to the general case. 
In the following figures we draw the graphs corresponding to figures 1-4 with 
k=0: 



- 2 2 -2 2 -2 2 

-2-1 1 2 

Graph corresponding to Figure 1. 



' In the infinite genus case this condition has to be weakened to the condition that 
the number is almost equal to (— 1)*2 in a sense, which depends on the class of 
solutions, which are considered (see [21]). 
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-2 2 -11 -2 2 



-2-1 1 



2 



Graph corresponding to Figure 2. 




Graph corresponding to Figure 3. 




Graph corresponding to Figure 4. 



Any neighbourhood of the point A = oo may be described by glueing together 
parts of finitely many sheets over some neighbourhoods of the point A = oo 
with all other complete sheets over A G CP^. Now v^— 1 In /i = arccos(A/2) 
defines a single-valued function on some of these neighbourhoods. Hence the 
glueing rules of all the parts of the sheets of some of these neighbourhoods 
must coincide with the corresponding glueing rules described by the standard 
graph. Now let us consider paths in the A-planc, which covers big circles in the 
A-plane moving clockwise around. They fit together to two paths with open 
ends. The corresponding succesion of those sheets, in which the path crosses 
the negative imaginary axis of the A-plane (and therefore the main part of 
the corresponding turn), results in two sequences of labels of our sheets: 

. . . , 2n + 2, 2n, 2n - 2, . . . and . . . , 2n - 1, 2n + 1, 2n + 3, . . . 



Hence for all graphs, which corresponds to periodic solutions of the non-linear 
Schrodinger equation, these paths should result in these two sequences of labels 
of our sheets. It is obvious that each part of these paths, which corresponds 
to one big circle moving clockwise around from the negative imaginary axis 
back to the negative imaginary axis in the A-plane results in a subgraph of 
the form indicated in figure 5, such that the following conditions are fulfilled: 
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(a) The real parts of the values of A at the branchpoints are ordered Re( Ai) < 
... <Rc{Ai). 

(b) Any sheet with label kj,j = 1, ... ,1 — 1 has no branchpoint, such that the 
real part of the corresponding value of A lies between Re(Aj) and Re(Aj_,_i). 

(c) The subgraph is maximal with respect to conditions (a) and (b). This 
means that the sheet with label ko has no branchpoint, such that the real 
part of the corresponding value of A is less than Re(Ai) and the sheet with 
label ki has no branchpoint, such that the real part of the corresponding 
value of A is bigger than Re(A;). 



. , A, . 

ko ki ki-i ki 

Figure 5. 

This gives the following condition on our graphs: 

(v) For each subgraph, which obeys condition (a), (b) and (c) the labels of 
the endpoints obey the equation ko — h — 2(— 1)^^° = 2(— 1)*^'. 

Now we can state the main theorem of this section: 

Theorem 3.2 The graphs described above are in one to one correspon- 
dence with the data, which corresponds to periodic solutions of the non-linear 
Schrddinger equation. 

Proof: Let us first prove that to all such graphs there exists some entire func- 
tion C — > C, A I— > A (A). In order to do this we divide the infinite number of 
cuts in three parts: To any part of the graphs of the form indicated in figure 
6. and figure 7. there corresponds the entire function 

C^C,A^ (-1)^2cos(VA) 



(-1)^2 (-1)^+^2 , (-1)^2 , (-1)^ +^2 

k k+1 k+2 k+3 

Figure 6. 



30 



... i-n 



(-1)^2 



(-1)^2 



k-2 



k-1 



k 



k+1 



Figure 7. 



Then due to conditions (ii) and (iv) the total covering space corresponding to 
some graph may be obtained by glueing in the middle of two such copies of CP^ 
a finite number of copies of CP^ along some cuts described by the glueing rule 
given by the graph. Due to conditions (i) and (iii) this total covering space 
is isomorphic to CP^. This shows that all these graphs define some infinite 
sheeted covering CP^ — > CP^ induced by some entire function A h- > A (A), 
where the parameter A is defined only up to some rational transformation, 
which leaves invariant the point oo. The transcendental curve defined by the 
equation /x^— /xA(A) + l = has infinitely many ordinary double points. In fact, 
all the branchpoints of order 1 of the covering map defined by these graphs, at 
which A is equal to ±2 correspond to such an ordinary double point. Hence 
the normalization of this transcendental curve is a compact Riemann surface. 
Due to condition (v) the function 1/ arccos(A/2) is a local parameter in the 
covering space near infinity. Hence there exists a unique global parameter A 
of the covering space, such that A = arccos(A/2) + 0(1/A) near infinity. This 
shows that the Riemann surface defined by the equation /x^ — /xA(A) + 1 = 
corresponds to periodic solutions of the non-linear Schrodinger equation. The 
other direction of the one to one correspondence is part of the definition of 
the graph. □ 

Corollary 3.3 The domain of the flow parameters {ti)i^z of the non-linear 
Schrodinger equation, defined in the first section may be chosen to be equal to 



More precisely, on this domain these flows have have only branchpoints of 
finite oder. 



In the proof we described implicitly how to determine the genus of the Rie- 
mann surface from the graph. In fact, to all branchpoints of the covering 

map induced by A there correspond two branchpoints of the covering map 
induced by fi from Y CP^. But those branchpoints, where A is equal to ±2 
correspond to ordinary double points and may be neglected, if they have no 
common sheet. 



(^i)iez £ all but a finite number ofti \ 



s are 



eqal to zero 



□ 
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The non-linear Schrodinger equation has two natural reality conditions. For 
both of them the function A (A) is assumed to be real, for all A G K. This 
condition implies that A(A) is equal to A(A). This gives the existence of an 
antilinear involution of the covering space of the covering map induced by 
A — > A(A) corresponding to the antilinear involution A — > A. Hence the cuts 
must be chosen to be invariant under this involution. For the starting point 
this may be achieved by turning the cuts: 

cut between the cut between the 

{2j — l)-th and the 2j-th and the 

2j-th sheet {2j + l)-th sheet 



There are more complicated situations, in which this antilinear involution 
permutes some of the sheets and the values of A are not real. In this case all 
but a finite number of cuts may be chosen to be invariant under the antilinear 
involution and therefore are part of the real axis A e R. But finitely many 
cuts are interchanged under this antilinear involution. We may choose them 
to be of the form 

A = Ao + and A = Aq - \^M.+ with ImAo > 0, respectively. 

The antilinear involution is then described by an involution of the sheets and 
the branchpoints, such that the value of A at two branchpoints, which are 
interchanged under this involution are the complex conjugate of each other. 
Let us give an example. In the following figure we describe the cuts by the 
pairs of the corresponding numbers of sheets. Again some number /c e Z is 
fixed: 



A = -2 



A = 2 



Figure 1'. 



(2j - 1, 2j) with j^k 



A = .5^^ 



(2j, 2j + 1) with j / fc, fc + 1 
and {2k -2,2k + 1) 



A = -2 




A = 2 



{2k- L,2fe + 1) 



Figure 8. 
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The involution of the sheets is given by the map 



Z — > Z, i I— > i for i 7^ A; — 1, A;; k — 1 i-^ k; and /c i— > A; — 1. 



Again we may combine this glueing rule to a graph: 



2 



-2 



2 



-2 



2 



2k-2 



21 



2k+2 



2k- 1 



Graph corresponding to Figure 8. 



Now we can describe the restrictions of the two reality conditions: The so 
called focusing case corresponds to the restriction that the values of A at the 
branchpoints are elements of [— oo, —2] U [2, oo]. In this case the links and ver- 
tices of the graph is always of the same form as the graph corresponding to 
figure 1. Then it is obvious that the coordinates (ti)iez are global single valued 
coordinates in the focusing case. The so called defocusing case corresponds to 
the restriction that the function A(A) takes values in [—2, 2] for all A e M. 
Hence the corresponding graph have to consist of one connected purely hori- 
zontal part combining all sheets and branchpoints, which are invariant under 
the involution describing the antilinear involution of the covering space of A. 
The corresponding values A have to be elements of [—2, 2] and have to be 
alternating with respect to the order. Beside this horizontal part the graph 
may have finitely many vertical strings of finite length, which are symmetric 
with respect to the horizontal part. This symmetry implements the antilinear 
involution. Hence the values of A above the horizontal part of the graph are 
assumed to have non-negative imaginary part and to be the complex conju- 
gate of the corresponding values below the horizontal part. In this case the 
coordinates (ti)iez are multivalued with complicated branchpoints. 

We want to remark that dX has a zero at the points where A^— 4 is equal to zero 
and d/i has a zeros at the branchpoints of the covering map induced by A. But 
all the ordinary double points may be neglected. Hence on the normalization 
of the algebraic curve defined by fi^ — fiA + 1 = the differentials dX and 
diJ, have a common zero, whenever the values of A at two branchpoints of a 
common sheet are equal to ±2. 

Finally let us mention that the differential equations describing the flows to 
the flow parameters (ti)iez can be solved numerically. We implemented these 
algebraic differential equations with the help of standard programs up to genus 
equal to 6. In particular, we checked the branchpoints of these flow parameters 
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numerically. With the help of such a program one can calculate quite quickly 
all the Riemann surfaces corresponding to the graphs described above. 
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